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ABSTRACT 

Of 104 AGN known to exhibit H2O maser emission, X-ray data that enable 
estimation of column densities, or lower limits, are available for 42. Contributing 
to this, we report analysis of new and archival X-ray data for 8 galaxies and 
collation of values for three more. Maser emission is indicative of large columns 
of cold gas, and in five of the eight new cases, maser spectra point toward origins 
in accretion disks viewed close-to edge-on (a.k.a. "disk- maser" systems). In 
these, we detect hard continuum and Fe Ka emission with equivalent widths on 
the order of 1 keV, which is consistent with Compton reflection, fluorescence by 
cold material, and obscuring columns > lO^^cm"^. Reviewing the full sample of 
42, 95% exhibit Nh > lO^^ cm'^ and 60% exhibit Nh > lO^^cm"^. Half of these 
are now recognized to be disk masers (up from 13); in this sub-sample, which is 
likely to be more homogeneous vis-a-vis the origin of maser emission, 76% exhibit 
Nh > lO^^cm"^. The probability of a common parent distribution of columns 
for disk- masers and other AGN masers is < 3%. Because ground-based surveys 
of AGN to detect new disk masers are relatively unbiased with respect to X- 
ray brightness and comparatively inexpensive, they may also be efficient guides 
for the sensitive pointed X-ray observations required to identify Compton-thick 
objects outside of shallow surveys. 

Subject headings: accretion disks — Galaxies: active — Galaxies: nuclei — 
masers — X-rays: galaxies 
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Introduction 



Ex tragalactic H 2 O ma s er emission a t 22 G Hz is known in 104 active galactic nuclei 
fAGNl Jzhang et all km^ . IXondratko I J2007h . ISraatz fc Gugliucci I J2008h . Braatz et al. 
2008, in preparation). In early analyses, the majority of the nuclei for which pointed X- 
ray observations exis t are Compton thi c k or nearly so, w ith columns ( Nh) ^ , 10^^ cm~^ 



flMadeiski et a/.ll2006l : IZhang et a/.ll2006l : iGreenhill Il2007h . IZhang et al\ mm consider all 
maser-host galaxies, where the emission is associated with either nuclear activity or star 
formation. In a sub-sample of 26 with apparent maser luminosity > 10 Lq (a.k.a. "mega- 



masers"), 85% exhibited Nh > 10^^ cm ^ and 50% exceeded 10^'' cm ^. iMadejski et al. 



(120061 ) focus on a probably more homogeneous sub-sample for which maser spectra suggest 
the emission specifically traces rotating, highly inclined disk structures close to the central 
engines. (We refer to these as "disk masers.") Seven of 9 nuclei with disk masers were 
Compton thick. 

It is plausible that AGN which host masers are characterized by high Nh in general, 
because large reservoirs of mole cular gas alon g the line of sight are required to produce 
substantial maser amplification (lElitzur Ill982l ). Among disk masers, prevalence of high Nh 
may be anticipated because the emission is believed to arise in inclined disk structures. 
We have analyzed new, archival, and published X-ray data for maser AGN and examine 



the st atistics of the updated sample, now about one third larger than that of IZhang et al. 



(I2OO6I ). Masers show promise as signposts of Compton-thick AGN that are readily detected 
from the ground. In §2, we describe new spectral analyses and estimation of Nh. We discuss 
selection effects in the full sample and inferences about maser AGN in §3. 



2. Data Analysis 

We observed seven maser-host AGN using the XMM-Newton observatory. In addition, 
we identified unpublished observations of five disk-maser hosts in the XMM and Chandra 
archives (12 targets total). All data were reduced using standard tasks in SAS and CIAO. 
XMM observations of NGC 6323 and NGC 6264 were lost due to flaring; NGC 3735 and 
NGC 2639 were lost due to high background. Table 1 lists the final sample. 

XMM EPIC PN spectra were extracted from 30" circles centered on the brightest pixels 
within 2" of the target coordinates (comparable to the 1.5" astrometric accuracy). Back- 
ground spectra were extracted from 80" circles, centered off-source, as far from the readout 
nodes as the targets. We did not use EPIC MOS data due to smaller effective area above 
2 keV. Chandra ACIS target spectra were extracted from 6" circles, while the backgrounds 
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were obtained from annuli centered on the targets, with radii of 15" and 30". The spectra 
were grouped with 10 counts bin~^ and analyzed using XSPECll. 

All 8 observations are characterized by relatively hard continua above 2keV. Impor- 
tantly, this is accompanied in five sources by Fe Ka emission with large equivalent widths 
(EWs), regardless of continuum model (Tab. 2; Fig.[T]). The strong lines probably arise from 
continuum reproces sing via fluorescence when the central engine is detected only in Comp- 



ton reflection (e.g., iLevenson et a/.ll2002l . 2006). Since signal to noise is relatively low for a 
majority of targets, we adopted a "minimum" model with a common set of fixed and free 
parameters. Because of degeneracies, results are best viewed as indicative. We assumed that 
the continuum is entirely described by reflection (a PEXRA V model) with underly ing (inci- 



dent) photon index fixed to 1.7. Larger values (F ~ 2; e.g.. IZdziarski et a/.lll996l ) observed 
toward type-1 AGN are also consistent. We modeled the Fe Ka lines as Gaussians with 
fixed 6.4 keV rest-frame energy and a = 20 eV intrinsic width. We also added a MEKAL 
component to match better emission below 2 keV. Fitted sub-solar abundances, when these 
arose, were probably not an accurate reflection of composition, but rather a consequence 
of a one-temperature model and omission of photoionization effects. Separate absorption 
columns were applied to MEKAL and PEXRAV plus line components. 

Spectral fltting with even the minimum model shows Fe Ka line emission in Mrk34, 
Mrkl419 and NGC6926, as well as marginally in NGC5793 (Tab. 2; Fig.[T]). We com- 
pute EWs with respect to the reflected continuum. In the case of NGC5793, we obtained 
EW=2.38 keV(x^/d.o.f=26.6/26) with a OkeV bound. Freeing the line energy, we obtained 
EW=2.46l|^^ keV and xVd.o.f.= 24.9/24, for line energy 6.32+°;°^ keV (rest). The shift in 
line energy may be due to coarse binning. Some elaboration on the minimum model was 
required for Mrk34, to accommodate excess near 7keV, which we flt with a second Gaus- 
sian at 7.05 keV (rest) with a = 20 eV, corresponding to Fe E Ws > IkeV for Mrk34 , 



Mrkl419, and NGC6926 are consistent with Nh > 1-2 xlO^'' cm'^ flLevenson et a/.l 120061 ). 
though quantitative interpretation depends on geometry. This and intrinsic limitations of 
the minimum model argue for inferred columns to be best regarded as indicative. No esti- 
mates are possible for ESO013-G012 and VV340A because neither exhibited clear Fe Ka 
emission (though both showed marginal excess at ~ 6 keV) . 
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Table 1. Observation Summary 



Galaxy 


Tel. 


Epoch^'^) 




Gounts 


Re 


NGC1320 


XMM 


06-08-06 


12260 


2778 


1 


NOG 1194 


XMM 


06-02-19 


12606 


1083 


2 


NGG6926 


XMM 


05-11-19 


8875 


149 


3 


Mrk 1419 


XMM 


05-10-30 


22298 


518 


4 


Mrk34 


XMM 


05-04-04 


8905 


682 


5 


NGG5793 


XMM 


05-07-21 


15686 


388 


6 


ESO013-G012 


XMM 


03-04-21 


21625 


167 




VV340A 


Ghandra 


06-12-17 


14862 


125 


7 



('^^Observing date in yymmdd ioirciaX,. 



•^''^ Duration of calibrated data used in modeling, in sec- 
onds. 

*^'^^Refererices fo r disk-maser identification in targets: (1) 
Kondratko I (120071). (2) (B raatz et ai, 2008, in preparat ion) 
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Greenhill et al 



Kondratko et al. 



Kondratko et al. 



2003h . f4l kondratko et a/.N2006h . (5) 



20061 
J2006h 



(6) 



Hariwara et al\ (119971 ). (7) 



Table 2. Summary of Spectral fitting 



Source 






kT 


Zq/Z ('')Fe Ka EW^-^) uh^'^) xVd-o.f. La^iokcV^'^ 




(lO^o cm-2) 




(keV) 


(keV) (10^4 cm-2) (ergs s-^) 



NGC 1320(^) 


65^^ 


1 42+0.55 


66+°-^^ 


[1.0] 


9 9n+0.44 
^•^U„Q 43 


> 1 


133/134 


7.5x10^^ 


NGC1194(^') 


99+?^ 


9 c:o+0.64 
^•^"-0.44 


91 +0-31 

"•^-'--0.26 


[1.0] 


9 1 1 +0.50 
^•-'--'--0.50 




69/74 


3.7x10^^ 


NGC 6926 


[7.43] 


[1.7] 


1 1 9+0.55 
-'-•-'-^-0.21 


[1.0] 


o OO+0.71 
O.Oo„3 21 


>i 


7.6/8 


5.8x10^° 


Mrk 1419 


[3.8] 


[1.7] 


0.76^°;- 


[0.05] 


1 7f^+0:89 
' "-0.88 


>i 


35/30 


3.2x10^° 


Mrk34('') 


[0.69] 


[1.7] 


0.791H^ 


[0.20] 




>i 


27/25 


9.8x10^1 


NGC 5793 


oq+12.8 


[1.7] 


0.481°:}^ 


[0.1] 






25/24 


1.8x10^^ 


ESO013-G012 


[7.39] 


[1.7] 


1 4Q+2.13 
-'-•^'J-0.51 


[0.1] 






5.7/11 


1.1x10^° 


VV340A 


90+28 


[1.7] 


U.UZ„g 25 


[1.0] 






5.1/6 


4.9x10^1 



Absorption due to Galactic and host-galaxy interstellar media. 

Values held fixed in fitting are enclosed in brackets. 

('^^ Computed with respect to the reflected continuum and for fixed 6.4 keV line energy (rest), except for the marginal 
detection in NGC 5793. 

Absorption within the AGN. Inferred indirectly from large Fe Ka EWs or directly from the absorbed power law 
component of NGC 1194. 

Apparent luminosity other than for NGC 1194, where intrinsic luminosity is 4.4 x lO^^ergss"^. 

^■^^Model includes a second MEKAL component with fixed temperature of 0.08 keV and a Gaussian component cor- 
responding to Fe K/? emission, fixed at 7.05 keV (rest) with intrinsic width a = 20 eV and fitted EW of 282+^^ eV. 



Normalization for the Fe K/3 line was constrained to be 0.11 x the normalization for the Fe Ka line ([George &: Fabian 



199lh . 

^^^As in (f) regarding secondary MEKAL and Fe K/3 components (EW=299l75 ^V). Additional Gaussian components 
representing blends of ionization states were fitted with observed energies of 1.75 ± 0.10 keV (Si Ka; EW=160;?^[]), 
2.49 ± 0.06 keV (S Ka; EW=390+^^°), and 3.58 ± 0.13 keV (Ca Ka; EW=380+^^°). The photon indices for incident 
reflected and transmitted components are set equal. 

('^^As in (f) regarding addition of an Fe K/3 component (EW=150'tgQ^ eV). 
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Comparatively high signal-to-noise ratios enabled fitting more complex models to data 
for NGC1194 and NGC1320. For NGC1320, we were able to fit a photon index, lA2toil, 
consistent with the value employed for other targets. The spectrum also shows excess consis- 
tent with Fe K/5 emission. In order to better represent the spectrum below 2 keV, we added 
a second MEKAL component with fixed 0.08 keV temperature. 

The spectrum of NGC1194 is still richer, exhibiting evidence above 5keV for direct 
transmission from the central engine, which we modeled as an absorbed power law with 



1.06 



+0.36 
-0.24 



X 10^^ cm ^ and F 



2 fi8+°-6^ 

^•00_Q 44. 



The column presumably arises in the 



reflector. About 70% of the estimated apparent source luminosity (2-10 keV) arises from 
transmission. Assuming that this component corresponds to emission from the central en- 
gine, the unabsorbed luminosity is 4.4x10^^ erg s~^, and the inferred Q/27T for the reflector 



i s on the order 0.1, which is consistent with t hat inferr ed for cold gas in, e.s 
( Colbert et aUbooi IPounds &: VaughanI [iooeh . Mrk 3 JPounds fc Pagelboosl ) 
JPone et a/.ll2003h . 



NGC 1068 
NGC 4945 



The spectrum of NGC 1194 also includes what appear to be Ka lines from Si, S, and 
Ca (Tab. 2). The profiles are broadened, which we hypothesize is d ue to photoionizatiq n 
and blending of emission from different states, as observed in IC 2560 (jMadejski et a/.ll2006l ). 
Because our primary intent is characterization of absorption columns, we only note in passing 
that treatment of photoio nization could be done in detail using an XSTAR model (e.g., 
Bautista fc Kallman |[200ll ). 



Discussion 



We have combined published estimates and limits on Nh with results of our analyses 
to obtain a sample of 42 AGN (out of 104 known maser hosts) . Th e sample includes: fl^ 



31 AGN that exhibit megam aser and weaker emission fjZhang et a/.ll2006l ). (2) two AGN 
hsted bv Izhang etall kood \ but without known columns - NGC 591, Mrkl, (3) five AGN 
treated here - Tab . 2, and (4) four n uclei in which maser emission has been recently dis- 
covered - NGC 17 kondratkoll2007h : NGC 3081, NGC 4253, and NGC 3783 (Braatz et ai, 
in preparation). Colu mns are taken from Zhan g et ai- 31 obje c ts, ex cept NGC 5256 (cf. 

Mrkl, NGC 17, 
NGC 3081; and 



Guainazzi et a/.M2005r): this work - 5 objects; iGuainazzi et a/.l (120051) 



and NGC 591 ; iReeves et al\ ( 120041 ) - NGC 3783; iMaiolino et all (119981 ) 



Pounds et al. J2OO3I ) - NGC 4253 In order of precedence, we adopt results from SAX, 



(Cappi et al. 


2006 


) to be more 


data alone {'. 


vomoto et al. 


2001) 
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Of the 42 AGN, 60% (25) are Compton t hick (Tab -El Fig. El). This fraction is consistent 

58% (18 of 31). Within 



with that obtained for the smaller sample of IZhang et al 



the current sample of 42, a greater fraction of recognized disk-maser galaxies are Compton 
thick, 76% (16 of 21). The distribution of columns for disk masers is significantly skewed 
and non-Gaussian (Fig.[2]). Substitution (over time) of Nh estimates for lower limits will 
exagerrate this. In contrast, the distribution for AGN masers without evidence of origins in 
edge-on disks includes a larger proportion of Compton-thin sources. A Kolmogorov-Smirnov 
(KS) test indicates a < 3.2% probability that the two distributions are drawn from a single 
parent. To account for lower limits in our test, we constructed a Monte Carlo simulation of 
1000 trials, sampled a uniformly distributed random variable of 10^^ and 10^^ cm~^ for each 
galaxy with a limit, and compiled a distribution of KS statistics. Boundaries of 10^^'^ and 
10^^ yielded probabilities of < 3% as well. 

Distinctions between Nh distributions for AGN that host acknowledged disk masers 
and those that do not may be still greater than estimated. Maser classifications depend on 
identification of spectroscopic markers: highly red and blue-shifted emission bracketing the 
systemic velocity. Disks may be oriented close to edge-on, resulting in large columns toward 
the central engine, while these spectroscopic markers may nonetheless be absent (e.g., high- 
velocity masers may be beamed away from the observer due to warping; cf. Miyoshi et al. 
1995). 



Zhang et al\ (120061 ) discussed an alternate sub-sample, selected based on apparent maser 
luminosity (i.e., megamasers). Within the current sample of 42 AGN with mas ers, we count 
58% (18 of 31) of recognized megamaser hosts are Compton thick (13 from IZhang et al. 



mm . plus NGC1194, NGC6926, Mrkl, Mrkl419, and Mrk34), matching the incidence 
among maser AGN broadly (Tab. [3]). We link the larger incidence of columns > lO^^cm"^ 
among disk-masers to a simple model geometry (edge-on orientation) and contrast this with 
selection based on maser luminosity, which is nearly always estimated with the convenient 
but inaccurate assumption of isotropic emission, without knowledge of beam angles and 
orientation to the line of sight. 

Columns > 10^^ cm~^ are more common in maser AGN than in a s ample of 22 hard X - 



ray (14-195 keV) selected, similarly nearby AGN observed with SWIFT ( iWinter et a/.ll2008l ). 
which includes no maser hosts (95% vs 50%; Tab. [3]). However, for Nh > 10^^ cm~^, incidence 
rates are similar (60% vs 50%). For masers in all classes of AG N, we estimate a statisti cally 
similar frequency of Compton-thick columns to that reported by lGuainazzi et al\ ( l2005l ) and 
Risahti et g/l fll999h among Seyfert 2 objects: 46% (22 of 48) and 65% (22 of 34), respectively. 
However, this is somewhat below the rate for masers associated with Seyfert 2 nuclei, 74% (20 
of 27). The difference is suggestive but comparison of column distributions of maser systems 



- 9 - 



Table 3. Column Measurements Among Known Maser AGN 



Sample > lO^^ cm'^ > lO^^ cm'^ No. 

Disk masers^'') 100% 76% 21 

All AGN masers 95% 60% 42 

Megamasers(^) 87% 58% 31 



^"^In addition to five sources in Tab. 2, we recognize 16 
other reported disk-maser AGN with accompanying esti- 
mates of Nh or limits: NGC591, NGC1068, NGC 1386, 
NGC2273, NGC 3079, NGC 3393, NGC 4051, NGC 4258, 
NGC 4388, NGC 4945, NGC 5728, IC2560, Mrkl, Mrkl210, 
Circinus, and 3C403. 

(^)We include 26 megamas e rs wit h estimates of Nh or lim- 
its as listed by IZhang et al\ (120061 ) . four objects in Tab. 2 
(NGC 1194, NGC 6926, Mrk34, Mrkl419), and Mrkl. 
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and local type-2 AGN is complicated by sample incompleteness in luminosity and distance, 
overlap between the X-ray and maser samples (e.g., 15/45 sources for Risaliti et al. 1999), 
and uncertainty in X-ray luminosities for obscured AGN (though aperture and extinction- 
corrected [OIII] luminosities could be used in place of X-ray estimates, e.g., Melendez et al. 
2008). 

Demonstrated correlation between incidence of maser emission and high Nh may suggest 
a new strategy for future efforts to identify obscured AGN. This is of particular interest if such 
objects are indeed substantial contributors to the Hard X-ray Background. Heavily obscured 
AGN typically have relatively low apparent luminosities and require long integrations (e.g., 
> 100 ksec to obtain > 10^ counts for cz = 10^ — lO'^kms"^). In contrast, maser emission 
in AGN as distant as cz ~ lO^kms"^ can be detected with integrations on the order of 
1 ksec, using a well-instrumented lOOm-dish antenna. We characterize the relative efficiency 
of maser and pointed X-ray observations, R, by the ratio of Compton-thick AGN discoverable 
in a fixed time by each means, for a sample of optically identified type-2 objects. We obtain 
i? ~ 3 from [10^/mMcT]/[10^Xc"r], where fm is the maser detection rate (~ 0.03; Kondratko 
et al. 2006), Mqt is the Compton thick fraction of maser AGN (~ 0.6; Tab. [3]), and Xct is 
the fraction of type-2 AGN found to be Compton thick in X-ray studies (~ 0.5; e.g., Risaliti 
et al. 1999). 

Deep, wide-field, hard X-ray surveys (e.g., NuSTAR) will extend synergy with maser 
studies, creating statistically significant X-ray and joint X-ray/radio samples. Nonetheless, 
there will be practical constraints on survey fields (5° x 5° for NuSTAR), and all-sky radio 
observations with well-defined completeness criteria could usefully guide pointed follow-up 
to enrich X-ray samples. 

Strengthened correlation between disk-maser emission and large obscurin g columns begs 
the question, is maser gas i n disks responsible for obscuration in type-2 objects (IGreenhill et al 



2003 



Madejski et a/.ll2006l )? On the one hand, obscuration by cold disk ga s at radii of 0.2-0.3 



pc is inferred in NGC 425 8 (IFruscione et al 



along other lines of sight (IGreenhill et al 



200 



20051 : iHerrnstein et al\ |2005| ) , and for Circinus 
On the other hand, time-variable Nh data 



i ndicate obscuration at radii of a fewxO.Ol pc toward NGC 1365, NGC 4388, and NGC 4151 
(lElvis et a/.ll2004l : iRisaliti et a/.l 120071 : iPuccetti et a/.l 120071 ). Direct physical association of 
maser and obscuring gas volumes still appears to be case-specific. Highly inclined structures 
that give rise to disk masers probably span a range of radii, at least extending inward with 
accretion fiows. Where these cross the line of sight to the central engine, and how extensively, 
likely depends on the detailed geometry of each system, requiring case study to establish. 



We thank A. Fruscione for help with data reduction and M. Elvis, A. Siemiginowska 
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Fig. 1.— XMM EPIC PN spectra of 8 AGN that host HgO masers. Each source was fitted 
with thermal emission model and a Compton reflection component (Tab. 2). All except 
VV340A and ESO013-G012 exhibit Fe Ka emission at 6.4keV (rest). Mrk34, NGC1320, 
and NGC1194 also show Fe K/5 at 7.05 keV (rest). NGC1194 exhibits evidence as well for 
(1) direct transmission of the central engine power law spectrum, obscured below ~ 5keV 
by a IMtoil X 10^^ cm ^ column and (2) Ka emission from Si, S, and Ca. 
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Fig. 2. — (left) Distribution of Nh in AGN that host masers and where X-ray data are 
available. Shaded bars indicate lower limits, {middle) Distribution for disk masers. Based 
on geometry, disk maser hosts are plausibly anticipated to constitute a more homogeneous 
sample with respect to Nh- (right) Distribution for maser AGN that are not identifiable as 
disk masers at this time. 



